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Abstract. In this work we take a critical look at the available data on the flyby 
anomaly and on the current limitations of attempts to develop an explanation. We aim 
to verify how conservative corrections to gravity could affect the hyperbolic trajectories 
of Earth flybys. We use ungravity-inspired potentials as a illustrative examples and 
show how the resulting orbital simulations differ from the observed anomaly. We also 
get constraints on the model parameters from the observed flyby velocity shifts. The 
conclusion is that no kind of conservative potential can be the cause of the flyby 
anomaly. 


PACS numbers: 04.50.Kd, 04.80.-y, 95.30.Sf 
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1 . Introduction 

Deep-space probes often use gravity assist manoeuvres to reduce the fuel requirements in 
order to reach their destinations. In the past few decades, some of those that used Earth 
gravity assists have shown what appears to be an unexpected shift in their hyperbolic 
excess velocities. This issue had already been the subject of discussion since the mid 
1990’s [H [2], but had maintained a relatively low prohle in the scholarly journals until 
a paper on the subject by Anderson et al. came out in 2008 [3]. This was at the height 
of the Pioneer anomaly controversy that had sparked a few years earlier BE], so the 
physics community was especially receptive to the discussion of spacecraft trajectory 
anomalies. The Pioneer anomaly was ultimately solved through conventional heat and 
radiative momentum transfer mechanisms mmn, with this solution being conhrmed 
by three independent efforts [HI El CO]. The same can also be stated about the Cassini 
anomaly observed while that spacecraft was cruising between Jupiter and Saturn, and 
whose cause is, most likely, also due to radiative momentum transfer m- However, 
unlike the Pioneer and the Cassini anomalies, this so-called flyby anomaly remains an 
open question ever since. 

This anomaly was detected in the residuals of the radiometric tracking data of 
several space probes performing Earth flybys. The trajectories inferred from the 
tracking data proved impossible to £t to a single hyperbolic arc. The pre-encounter 
and post-encounter data had to be £t by separate hyperbolic trajectories that displayed 
a discrepancy in the hyperbolic excess velocities. Since the two hyperbolic arcs appear 
to be good hts for their respective data sets, it is assumed that this velocity shift is 
localised near the perigee, where tracking through the Deep Space Network (DSN) is 
not available for approximately four hours [3]. This fact has, indeed, led to a proposed 
experimental setup to improve this time window [12]. The spatial resolution of the 
available reconstructions, resulting from the 10 s interval tracking, does not allow for 
an accurate characterisation of the effect, so that no corresponding spatial or temporal 
prohle of the acceleration exists. The only available data to support the analysis is the 
shift in the excess velocity (and correspondingly, in kinetic energy) between the inbound 
and outbound trajectories. 

The hyby anomaly has been observed in the Galileo, NEAR, Rosetta and Cassini 
hybys. Earth hybys between 1990 and 2005 are listed in Table [H based on data from 
Ref. [3], with the respective perigee altitudes (hp) and velocities (up), hyperbolic excess 
velocities (uoo) and the anomalous shift in the excess velocity (Auoo)- 

It should be noted that the information on the hybys where the anomaly was 
allegedly detected is not entirely consistent across the diherent sources. Aside from 
minor disparities in the values of some parameters, there are a few important diherences 
that deserve some discussion. 

Regarding the second Galileo hyby. Ref. [3] explains that the measured value for 
the change in excess velocity (Auoo) was actually —8 mm/s and the —4.6 mm/s hgure 
in Table [T] is obtained after subtracting an estimate of the atmospheric drag. However, 
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Table 1. List of Earth flybys between 1990 and 2005, based on data from Ref. [3], 
where hp and Vp denote the altitude and velocity at the perigee, Voo is the hyperbolic 
excess velocity and Avoo is the measured change in excess velocity. 


Date 

Mission 

hp 

(km) 

Vp 

(km/s) 

(km/s) 

(mm/s) 

08/12/1990 

Galileo 

960 

13.740 

8.949 

3.92 

08/12/1992 

Galileo 

303 

14.080 

8.877 

-4.6 

23/01/1998 

NEAR 

539 

12.739 

6.851 

13.46 

18/08/1999 

Gassini 

1175 

19.026 

16.01 

-2 

04/03/2005 

Rosetta 

1956 

10.517 

3.863 

1.80 

02/08/2005 

MESSENGER 

2347 

10.389 

4.056 

0.02 


for the same flyby, two other sources claim that, due to the low altitude, atmospheric 
drag would mask any anomalous velocity change mm- The reasons for absence of any 
value for the August 1999 Cassini flyby in Ref. |T3] are also uncertain, as no explanation 
is offered for this omission. 

When one compared the tables presented in Refs. [3j and [13], one is not merely 
facing slight differences in values. The two sets of data in these two papers could actually 
imply a whole different picture on the phenomenon, potentially prompting different 
assumptions. While in the hrst the effect appears to be bidirectional, accelerating the 
spacecraft on some occasions and decelerating in others, in the second reference all 
instances that have negative velocity changes have been removed and the effect appears 
to be consistent with an energy increase. 

The study of more flybys and closure of the four hour gap in DSN tracking near the 
perigee would be essential to shed some more light onto this phenomenon. There have 
been a few recent flybys, listed in Table [2], however their results are yet to be presented 
in scholarly literature. 

Table 2. Recent Earth flybys for which no data is yet available in the literature. 


Date 

Mission 

13/11/2007 

Rosetta 

13/11/2009 

Rosetta 

09/10/2013 

Juno 

03/12/2015 

Hayabusa 2 


Some attempts have been made to hnd explanations for the flyby anomaly, or 
empirical descriptions like the one put forward in Ref. [5]. A very detailed of the 
two Galileo (1990 and 1992) and the NEAR (1998) gravity assists considered effects 
of Earth oblateness, other Solar System bodies, relativistic corrections, atmospheric 
drag. Earth albedo and infrared emissions, ocean tides and solar pressure [1]. This 
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analysis was extended to other possible error sonrces snch as the atmosphere, ocean 
tides, solid tides, spacecraft charging, magnetic moments, solar wind and spin-rotation 
conpling m. None of these efforts was able to hnd a snitable explanation. Specifhcally, 
the effect of Earth oblateness could yield an acceleration with a compatible order of 
magnitude, although all attempts yielded unreasonable solutions, unable to account for 
all flybys |1]. A study of the thermal effects in the hrst Rosetta flyby has also been 
preformed, concluding that they cannot be responsible for the reported flyby anomaly 
[9]. A discussion about some exotic explanations can be found in Refs. [ElUS], including 
interaction with dark matter [16], modihed inertia Id or modihed particle dynamics 
H. There are other efforts underway, for instance at the lAU, looking into the effects 
of the dehnition of coordinate frames on hyby trajectories |18j . 

At this point, it must be acknowledged that the hyby anomaly is still poorly 
characterised and any consistent treatment is, at least, very difficult. These limitations 
clearly raise the need for alternative approaches. 

In this paper, we hrst perform an order of magnitude analysis on some mechanical 
ehects that could putatively account for the hyby anomaly. This allows us to acquire 
some sensitivity on the nature of the phenomenon. We then use an ungravity-inspired 
potential as an example of a conservative modihcation to gravity and study its ehects on 
the Earth hyby trajectories. This allows us to reach general conclusions on the ehects of 
conservative potentials on the suspected anomalous hyperbolic trajectories. The main 
idea to retain is that conservative forces cannot explain the hyby anomaly. 

2. Newtonian Dynamics 

The data available about the hyby anomaly suggests there may be an energy (and linear 
momentum) shift highly localised at the perigee. It is sensible to look for conventional 
means for that shift to take place. This procedure also allows one to gain some sensibility 
on the hgures involved. 

We shall thus begin by examining some conventional Newtonian mechanisms, if 
nothing else, to acquire some sensitivity on the hgures of merit. This analysis shall be 
taken in terms of orders of magnitude. 

The most obvious mechanism to induce a change in orbital energy is the separation 
of a certain amount of mass with a certain speed relative to the spacecraft, in ehect, 
a thrust. In this discussion, we shall not be concerned with the specihc origin of this 
mass loss, but instead to attempt to estimate an order of magnitude for the speed and 
amount of mass that would have to be involved in order to generate a shift in energy 
similar to the one reported in the hyby anomaly. 

The hrst Galileo hyby of December 8, 1990 provides a good benchmark for an 
order of magnitude analysis. It had a perigee altitude of approximately 10^ km with a 
maximum velocity of around 14 km/s. Its hyperbolic excess velocity was Uoo ~ 9 km/s, 
to which corresponds a specihc orbital energy (energy per unit mass) T ~ 4 x 10^ J/kg. 

The frequency shift implied by the residuals for this hyby could be explained by an 
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increase in hyperbolic excess velocity Awoo ~ 4 nini/s or, equivalently, a shift in specihc 
energy A£ k, 35 J/kg. This corresponds to a An 2.6 mm/s at the perigee. 

Assuming that this An is the result of some kind of thrust, we can, from basic 
mechanics, through conservation of linear momentum, obtain a relation between the 
mass and speed of the propellant. It is not relevant for this analysis to speculate if this 
propellant is a solid fragment or a gas. We consider the most favorable case of thrust 
along the flight direction yielding the largest possible effect an thus obtaining an upper 
bound of this effect. Also, for simplicity, we assume the thrust to be instantaneous, which 
seems reasonable given the highly localised nature of the anomaly we are searching for. 

The results we hnd from this analysis show that any fragment or amount of mass 
small enough to be accidentally lost without being detected, would have to have a 
speed relative to the spacecraft itself of, at least, the same order of magnitude as the 
spacecraft’s velocity relative to Earth. Indeed, for a fragment or leak small enough for 
an undetected loss, let us say of around 10“^ of the mass of the spacecraft or 0.25 g, 
the required relative speed would be of the order of 10^ m/s for an ejection along the 
direction of the motion, the most favourable angle. 

Even though there is no known internal mechanism for such thing to happen without 
being detected, one could consider the reverse process. If the spacecraft were to collide 
with a piece of debris or a meteoroid, relative speeds of this order are, indeed, possible. 
If the collisions were to take place at the right angles, this could be a possibility to obtain 
the kinds of velocity shifts obtained, although the spacecraft could sustain significant 
damage. However, a recently published survey of damage sustained by the Space Shuttle 
from micrometeoroids and orbital debris has found that the most frequent damaging 
fragments have masses ranging from 10“® g to 10“^ g [19], which is at least two orders 
of magnitude below what we used in our estimate, for comparable speeds. Meteoroids 
with the required mass must be rare enough to be ruled out as an explanation for the 
flyby anomaly. 

Another speculation to achieve the reported energy shift is a change in rotational 
kinetic energy resulting from a transfer of translational kinetic energy, with the 
consequent coupling between linear and angular momentum. Again, we do not attempt 
to propose a specihc mechanism for this, but only to assess the viability of such a 
putative process through an estimate of the involved orders of magnitude of the energy 
exchange. 

From the approximate dimensions of these typical spacecraft {e.g. Galileo is a 
~ 300 kg spacecraft with an approximately cylindrical shape with a ~ 1 m radius, a 
~ 3 m height), one can estimate the moment of inertia of the spacecraft and compute 
the change in angular velocity due to the translational energy change. A typical hgure 
of merit is around 10 rad/s, which is clearly too large to go undetected. Thus, one 
can conclude that an explanation for the flyby anomaly from any kind of rotational 
mechanism is unfeasible. 

This analysis, even though straightforward and approximate, alows for acquiring 
some sensitivity on the hgures involved in the physical processes, from a classical 
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mechanical standpoint. Hence, we verify that, despite being a small effect from the 
orbital analysis point of view, the flyby anomaly involves a shift in velocity and energy 
that is too large to be caused by any kind of subtle hypothetical mechanism. 

3. An example of a conservative modification to gravity: 

Ungravity-inspired potential 

3.1. Unparticles 

Recently, the possibility of the existence of new physics above the TeV scale has been 
considered through the introduction of unparticles [20l [21]. In this scheme one admits 
a hidden sector with a nontrivial infrared hxed point Au, below which scale invariance 
is explicit. In the ultraviolet (UV) regime, at energies above Ay, the hidden sector 
operator C>uv of dimension duv couples to the standard model (SM) helds through an 
operator Osm of dimension n via nonrenormalizable interactions , 

where Mu is the mass of the heavy exchanged particle. Below Au, the hidden sector 
becomes scale invariant and the operator Ouv uiutates into an unparticle operator 0\] 
with noninteger scaling dimension d^. The coupling of held operators can be generically 
written as 

A C^UV—<iu 

^rfuv+n-4^uC)sM (1) 

The operator 0\] could be a scalar, a vector, a tensor or even a spinor. 

The exchange of unparticles gives rise to long range forces which deviate from 
the inverse-square law (ISL) for massless particles due to the anomalous scaling of the 
unparticle propagator. For example, the exchange of scalar (pseudo-scalar) unparticles 
can give rise to spin-dependent long range forces, as pointed out in Ref. |22] . Coupling 
between unparticles and vector or axial-vector currents have been investigated in Ref. 
|23] . In Ref. [23] the coupling between unparticles and the energy-momentum tensor 
was studied. An analysis of some of the phenomenological implications can be found in 
Refs. [25l|26l|27]. 


3.2. Ungravity 


If Ou is a rank-2 tensor, it can couple with the stress-energy tensor and lead to 
a modihcation to Newtonian gravity. Taking the gravitational coupling of the tensor 
unparticle to to be of the form 

( 2 ) 


MA' 




where M* = Au{Mu/AuY^^, it can be shown, in the non-relativistic limit and for 
du 7^ 1, fhat the effective gravitational potential with the unparticle exchange has the 
form [23] 

2(iu-2' 

/ \ 

1 + 


V{r) = 


GnM 


r 


r 


(3) 
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where is the Newtonian gravitational constant and Rq is a characteristic length 
scale. 

One can also obtain a force from the coupling of a vector unparticle [23]. The 
potential for the coupling between a vector unparticle and a baryonic (or leptonic) 
current of the form 


A 


- 1 

u 


is given, when combined with the gravitational potential, by 


(4) 


V{r) 


GnM 

r 



(5) 


In this paper, we use Eqs. (jS]) and (jS]) from a phenomenological model standpoint, 
and we shall refer to them, respectively, as “tensorial” and “vectorial” ungravity-inspired 
potentials. In this context, we attempt to constrain the values of Rq and d^ from the 
velocity and energy shifts observed in the Earth flybys discussed in Section [1] This is 
analysed from two different perspectives. First, we look at the way the hyperbolic 
trajectories of these spacecraft would be altered by the existence of an ungravity- 
inspired potential and, particularly, what is the temporal signature of the radial velocity 
perturbation. In a second analysis, we attempt to obtain the range of values for the 
potential paramenters Rq and du that leads to a velocity anomaly that hts the order of 
magnitude of the observed instances. 


4. Results and Discussion 

The existence of a potential of the kind shown in Section [3] would have observable effects 
on the trajectories of objects under gravity. Specihcally, the hyperbolic trajectories of 
the Earth flybys under analysis in this paper would show a deviation from their predicted 
trajectories. 

In this Section, we aim to characterise the perturbation induced by an ungravity- 
inspired potential of the form shown in Eqs. ([3|) and ([5|). In order to achieve that, we 
have performed numerical simulations of the orbits reproducing the conditions of the 
flybys and repeated them with the addition of the new potential. We can then compare 
the two trajectories in the search for the observational signature of the perturbation and 
the bounds on the parameters established by the flyby data. 

4.I. Temporal Signature 

The anomalous velocity shift that was discovered in some Earth flybys is inferred from 
the inability to £t the inbound and outbound arms of the trajectory to a single hyperbolic 
segment. It shows up in the Doppler data residuals as a sudden step in frequency at 
the perigee. 

The Doppler frequency shift is a function of the radial velocity of the probe. For this 
reason, the easiest way to look for a Doppler frequency shift from a trajectory simulation 
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dr(t) 

A- 

dt 



Figure 1. Difference between the radial velocity of the hypothetical trajectory with 
a tensorial ungravity-inspired potential and the Newtonian trajectory as a function of 
time. Both trajectories are simulated for the same initial conditions. The tensorial 
ungravity-inpired potential uses Rq = Rq and du = 1 -I- (5 x 10“^), 

is to extract the radial velocity perturbation resulting from the new potential and look 
at its time signature. This is a similar methodology to the one used in Ref. |28] . 

At this stage, we will not look at the values of the parameters, though some 
preliminary simulations have shown us that for values of Rq near the Earth radius, 
R®, the order of magnitude of du must be within 1 ± 10“® in order to produce a velocity 
shift that is remotely compatible with the ones reported in the flyby anomaly. 

When we plot the perturbation in radial velocity relative to the Newtonian 
trajectory as a function of time, the most striking feature is a large spike localised 
between approximately 10 minutes before and after the perigee, as can be seen in Figs. 
[U and |2] . Besides that, there is a smaller deviation in the inbound and outbound 
trajectories, but that converges to the Newtonian value as we get further from the 
perigee. Figs. [1] and |2] also show the difference between applying a “tensorial” or 
“vectorial” unparticle potential, which basically flips the graphic in the vertical direction. 

From the observation of these results, it is already obvious that an ungravity- 
inspired potentials cannot produce a perturbation that would yield the kind of temporal 
signature that is observed in the Doppler residuals of the Earth flybys. Still, one might 
be tempted to reason that, if the large spike around the perigee would somehow not 
appear in observational data due to lack of resolution, what remains is a step in radial 
velocity that somewhat resembles what is observed in the residuals [3]. One should, 
thus, look closely at Figs. [H and [2] and observe, once again, that the velocity shift 
quickly converges to zero on both the descending and ascending trajectories as one gets 
further away from the perigee, a behaviour consistent with energy conservation, unlike 
the flat behaviour that is observed in the residuals that implies a non-conservative effect. 





Hyperbolic orbits of Earth flybys and effects of ungravity-inspired conservative potentialslO 


dr(t) 

A- 

dt 



Figure 2. Same as Fig. [U for a vectorial ungravity-inspired potential, with Rq = 
and du = 1 -I- (5 X 10”"^). 

Furthermore, the way in which these time signatures were obtained means that they are 
not directly comparable to the Doppler residuals. 

We can stress this point further by attempting to reproduce the conditions in which 
those residuals were produced. To do that, we now compare two Newtonian hyperbolic 
arms obtained separately from the inbound and outbound parts of an underlying 
hypothetical trajectory with the ungravity-inspired potential. 

In this case the results show a similar spike in the radial velocity shift near the 
perigee, although the overall time signature is now symmetric around the perigee. This 
last feature is due to the fact that the perturbation to the potential is still an inverse 
power law of the radius. This means that the integration of the perturbation along 
any time interval centred at the perigee will always lead to a symmetric time signature 
converging to zero at the extremities. 

From this last observation, it is logical to conclude that, since any kind of potential 
based on position dependent deviation to the inverse square law will be unable to explain 
the flyby anomaly. One can reinforce this argument by recalling that a potential that 
preserves energy conservation can only depend on the coordinates. Since the velocity 
shift in the flyby anomaly implies a corresponding orbital energy shift, it cannot be 
explained by any kind of conservative force. 

4-2. Constraints on Rq and d^ 

Although it is now clear that the temporal signature of ungravity-inspired potentials 
are not compatible with the one reported in the flyby anomaly, it is still interesting to 
attempt to use the velocity shifts involved in these flybys to constrain the ungravity 
parameters Rq and d^. 

To bound these parameters, we had to perform a series of numerical simulations 
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Figure 3. Difference between the radial velocity of the outbound and inbound 
Newtonian trajectory arms computed from the underlying hypothetical trajectory 
perturbed with a tensorial ungravity-inspired potential as a function of time. This 
process aims to reproduce the way in which residuals are obtained in Ref. [3]. The 
tensorial ungravity-inpired potential uses Rq = Rq and du = 1 -|- (5 x and the 
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Figure 4. Same as Fig. [31 for a vectorial ungravity-inspired potential, with Rq = Ri 
and du = 1 + (5 X 10“”^). 









Hyperbolic orbits of Earth flybys and effects of ungravity-inspired conservative potentialsl2 



Figure 5. Contour plot of the maximum radial velocity shift, in m/s, of the trajectory 
perturbed by the ungravity potential relative to the Newtonian trajectory, as a function 
of i?G and d-a- The plot is similar for tensorial and vectorial potentials. 


where the maximum radial velocity shift caused by the ungravity potential is measured. 
We then plot the maximum velocity shift as a function of Rq and d^- Due to the 
characteristic perigee altitudes involved in these flybys, we look for characteristic length 
scales near the Earth radius, i.e., Rq ~ R^. A preliminary set of simulations also 
indicated that the values of du should be within 10“® of unity. 

The result is depicted in Fig. [5]where the dependence of the maximum perturbation 
shift on the ungravity parameters is shown. In this case, we tested values of d^ below and 
above unity, though always very close to it, as discussed above. The results are similar 
for tensorial and vectorial ungravity, since, as we have also discussed, the difference 
between these two only causes the vertical flipping of the time signatures, which is the 
same effect that the change from du > 1 to du ^ 1 has. 

What these results show is that the perturbation grows as du is further away from 
1, as would be expected. The same can be stated as Rq gets smaller. More importantly, 
we can see that the range of values that is compatible with a velocity shift of the order 
of 10 mm/s or less is a narrow strip around du = 1, that gets wider as Rq grows. Still, 
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a much larger Rq would start to affect other orbits at larger scales and would not yield 
a localised effect like the flyby anomaly. 

Overall, we stress once again that, not only do the flyby anomaly measurements 
place very stringent limits on but also that this parametric analysis does not take into 
account the above discussion about the time signature of the effect, where we pointed 
out that no conservative effect can alone explain the flyby anomaly. 

As a hnal remark, we could point out that, if the effect of a conservative force to be 
added to the Newtonian component changes the perigee altitude, then the importance 
of dissipative forces such as atmospheric drag might be expected to change as well, as 
density has an exponential dependence with altitude. However, for the typical altitudes 
involved, the density is small enough so that the effect of atmospheric drag on the 
velocity is already much smaller that the anomaly. Clearly, this is not relevant for the 
cases depicted in Tabled] but the fact remains that the only possibility for a conservative 
force to be relevant in the flyby anomaly is through a coupling with some kind of non¬ 
conservative force. 

5. Conclusions 

The flyby anomaly has been a puzzle for quite some time and, so far, no credible 
explanation has emerged. This paper does not provide a hnal answer either, but 
eliminates some possible causes and raises attention for the time signature test that 
any possible explanation of the hyby anomaly must pass. 

Given the discussion we have made in Section dl a signihcant effort must still 
be made in what concerns the characterisation of the phenomenon. It would be 
extremely important to have results on the trajectory analysis for the most recent hybys. 
Otherwise, the six observational instances on which this whole discussion is based, with 
no known independent analyses and with all the caveats that we have discussed, may not 
be sufficient to unequivocally establish the features that are attributed to this anomaly, 
or indeed its existence. The proposal of Ref. [12] can provide a relevant contribution in 
this direction. 

The discussion around these ungravity-inspired potentials sets out another example 
of an attempt to tackle this anomaly based on a modihed gravity model. It gives 
a characterisation of the time signature that this kind of modihcation to Newtonian 
gravity yields. For the parameter range tested in this paper, namely, for length scales 
similar to the Earth radius, we obtain a highly localised effect near the perigee. Still, the 
effect is symmetric around the perigee, with the perturbation converging back to zero as 
the object goes further from the Earth flyby, as would be expected for any conservative 
force. We also constrained the length scale and non-integer dimension of the ungravity 
potential with bounds set by the flybys where the anomaly was detected, leading to a 
very narrow range of allowed values for this last parameter. 

In the end, we must conclude that an potential inspired by unparticle-like 
corrections, such as the ones used in this work, cannot explain the flyby anomaly. 
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Indeed, the observations made abont its time signatnre remain valid for any inverse 
power law of distance which necessarily conserve orbital energy. The flyby anomaly has 
an inherently non-conservative nature. 

The search for a solution for the flyby anomaly has been akin to the work of a 
detective, and has yet to provide any palpable results or, indeed, signihcant advances 
in the knowledge of the anomaly itself. Until more observational data is available, we 
predict that the situation is unlikely to change. 
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